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(U)  The  probability  density  function  (PDF)  for 
linear  glint  has  been  calculated  for  the  situation  in 
which  the  tracking  angles  derived  from  the  co-  and 
cross-polarized  echoes  are  combined  to  produce  an 
averaged  tracking  angle.  A  target  consisting  of  two 
scattering  centers,  each  producing  an  echo  signal  with 
a  Rayleigh  distribution,  was  assumed,  and  two  methods 
of  combining  the  tracking  angles  were  investigated. 

The  PDFs  were  generated  using  a  Monte  Carlo  simulation. 
The  results  show  that  combining  the  two  tracking 
angles  in  this  manner  can  reduce  the  fraction  of  the 
glint  PDF  falling  outside  the  bounds  of  the  target 
from  0.26  (the  fraction  when  either  echo  is  used  by 
itself)  to  0.11,  if  the  co-  and  cross-polarized  echoes 
are  not  correlated.  If  the  echoes  are  correlated,  the 
improvement  varies  linearly  with  the  correlation 
coefficient  between  these  two  values.  Combining  using 
weighting  of  the  derived  tracking  angle  by  the  radar 
cross  section  in  each  echo  produces  a  better  improve¬ 
ment  than  simple  averaging  of  the  angles...  - 
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INTRODUCTION 


Tracking  radars  follow  a  target  by  keeping  the  antenna  oriented 
normal  to  the  echo  wave  front.  For  relatively  distant  targets  or  for 
simple  targets  with  one  dominant  scattering  center,  the  wave  fronts  are 
spherical,  with  the  center  of  the  sphere  located  at  the  target,  so  that 
the  antenna  indeed  points  toward  the  physical  target  center.  However, 
at  close  ranges  and  for  more  complex  distributed  targets,  interference 
between  reflections  from  two  (or  more)  strong  scattering  centers  on  the 
target  produces  distortions  in  the  phase  fronts.  The  phase  front  dis¬ 
tortions  can  easily  be  severe  enough  that  the  radar  points  many  degrees 
away  from  the  true  center  of  the  target  in  an  attempt  to  maintain 
normality  to  the  local  phase  front. 

In  some  tracking  applications,  such  angle  scintillation  can  be  a 
severe  problem,  and  various  mitigation  techniques  have  been  studied. 

Since  the  interference  is  strongly  dependent  on  frequency  and  the  spatial 
location  of  the  receiving  antenna,  frequency  and  spatial  diversity  have 
been  proposed  and  investigated  by  analytical,  simulation,  and  measurement 
methods.  Our  interest  has  focused  on  the  use  of  polarization  diversity 
as  a  mitigation  technique,  under  the  assumption  that  the  phase  and 
amplitude  of  the  cross-polarized  return  is  significantly  decorrelated 
from  the  co-polarized  return.  Poelman1  and  Varshavchuk  and  Kobak,2 
among  others,  have  shown  that  the  correlation  between  orthogonally 
polarized  backscattered  components  of  an  extended  target  can  vary  from 
0  to  1,  depending  on  details  of  the  target,  incident  wave  polarization, 
etc.  In  this  report,  we  incorporate  the  correlation  as  a  parameter. 
Polarization  diversity  may  be  more  readily  implemented  in  some  appli¬ 
cations  than  either  frequency  or  spatial  diversity. 

A  first  step  in  assessing  the  utility  of  polarization  diversity  is 
to  compute  the  probability  density  function  (PDF)  for  various  methods  of 
combining  the  co-  and  cross-polarized  returns,  as  a  function  of  the 

JA.  J.  Poelman.  "Cross  Correlation  of  Orthogonally  Polarized 
Backscatter  Components,"  IEEE  Trans.  Aerospace  anti  Electron.  Sue  terns, 

Vol.  AES-12  (1976),  pp.  674-681. 

2M.  L.  Varshavchuk  and  V.  0.  Kobak.  "Cross  Correlation  of 
Orthogonally  Polarized  Components  of  Electromagnetic  Field  Scattered  by 
an  Extended  Object,"  Radio  Eng.  Electron.  Phys.  (USSR),  Vol.  16  (1971), 

pp.  201-205. 
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correlation  between  the  two  returns.  The  purpose  of  this  report  is  to 
present  the  results  of  such  a  calculation.  We  assume  a  relatively 
simple  target  with  two  isotropic  scattering  centers,  each  with  a 
Rayleigh-distributed  radar  cross  section  (RCS) .  We  investigate  two 
combining  techniques:  simple  averaging,  and  averaging  in  which  each 
derived  tracking  angle  is  weighted  by  the  RCS. 

Because  of  analytical  difficulties,  a  closed  expression  cannot  be 
derived  for  the  PDF  of  the  combined  returns,  and  Monte  Carlo  techniques 
must  be  used.  A  Monte  Carlo  technique  also  has  the  advantage  that  the 
correlation  between  the  co-  and  cross-polarization  can  be  incorporated 
in  a  straightforward  manner. 

The  results  presented  indicate  that  diversity  combining  of  the  co- 
and  cross-polarization  returns  can  reduce  the  probability  that  the 
tracking  angle  falls  outside  the  boundaries  of  the  two-scatterer  target 
from  0.26  to  0.11.  The  reduction  depends  on  the  correlation  of  the  co- 
and  cross-polarization  echoes,  and  the  reduction  is  greatest  for  the 
RCS-weighted  combining.  The  basic  conclusions  are  not  expected  to  be 
changed  by  more  complicated  targets  with  more  scattering  centers. 


THEORY 


The  geometry  of  the  target  and  tracking  radar  is  shown  in  Figure  1. 
The  target  consists  of  two  independent  scattering  centers  whose 
scattered  electric  field  amplitude  distributions  are  described  by 
Rayleigh  PDFs: 


f(E 


)exP(-Ei>j/2a^j) 


(1) 


where  i  =  1,2  is  an  index  denoting  the  scattering  center;  j  =  co,cr  is 
an  index  labeling  the  co-  or  cross-polarized  scattered  field,  respec¬ 
tively;  and  o  is  the  standard  deviation  of  the  Rayleigh  PDF.  As 

1  9  J 

derived  by  the  tracking  loop  in  either  a  conical-scan  or  a  monopulse 
radar,  the  apparent  radar  center  of  this  two-scatterer  target  is 
displaced  from  the  true  geometrical  center  by  a  distance  R  because  of 
the  phase  front  distortion.  For  two  scatterers  separated  by  a  distance 
D,  the  value  of  R  derived  from  the  co-polarized  return  is  given  by 

R  -  (i)(l  -  Z2  )(1  +  2Z  cosiJj  +  Z2  r1  ,  (2) 
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FIGURE  1.  Target  and  Seeker  Geometry. 

where  ZCQ  »  ElfCO/E2,Co  and  ^co  is  the  Phase  difference  between  Ej >c0 
and  E2  co.  The  phase  difference  can  arise  both  from  a  difference  in 
path  length  from  each  scattering  center  to  the  receiver  or  from  phase 
shifts  caused  by  the  scattering  centers.  A  similar  expression  holds 
for  Rcr  for  the  cross-polarized  return  in  terms  of  Zcr  and  if»cr.  This 
equation  was  first  derived  by  Locke3  and  later  by  other  approaches  by 
Ostrovityanov,1*  Dunn  and  Howard,5  and  Dunlop.6  As  shown  in  Reference  6 
the  equation  requires  an  automatic  gain  control  (AGC)  loop  with  a  time 
constant  low  enough  that  glint  fluctuations  are  tracked  instantaneously 
(but  with  a  time  constant  slower  than  the  scan  rate  for  a  conical-scan 
system) . 


3A.  S.  Locke.  Guidance.  New  York,  Van  Nostrand  Reinhold,  1955. 

4R.  V.  Ostrovityanov.  "Angular  Noise,"  Radiotekh.  Eleckron.  , 

Vol .  4  (1966),  pp.  507-515. 

5J.  H.  Dunn  and  D.  D.  Howard.  "Radar  Target  Amplitude,  Angle, 
and  Doppler  Scintillation  from  Analysis  of  the  Echo  Signal  Propagating 
in  Space,"  IEEE  Trans.  Microwave  Theory  Tech.,  Vol.  MTT-16  (1968), 
pp.  715-728. 

6A.  J.  Dunlop.  "AGC  Response  and  Target  Glint,"  Proc.  IEE , 

Vol.  128,  Pt.  F,  No.  2  (1981),  pp.  83-90. 
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Equation  2  is  plotted  in  Figure  2.  Phase  differences  near 
180  degrees  and  amplitude  ratios  near  1  are  observed  to  produce  an 
apparent  radar  center  that  falls  outside  the  boundaries  of  the  target. 


0 


1  FIGURE  2.  Plot  of  Equation  2.  The  apparent  radar  center  is  given  in 

terms  of  the  scatterer  separation.  The  positions  of  the  scatterers  are 
|  indicated  by  the  dotted  lines. 


The  PDF  for  Rco  for  two  scattering  centers,  each  with  Rayleigh  PDFs 
and  a  phase  difference  i|>c0  uniformly  distributed  over  0  to  180  degrees, 
has  been  derived.4  The  PDF  is  given  by 


2y 


f(R  ) 

CO 


CO 


r. 


Y*  [<1  + 
CO 


2R  o  o 

IT*  +  4(1 


2RcOy2  3/2 
D  '  J 


(3) 


where  yc0  =  oi  co/a2  co*  For  the  Present  calculation,  the  PDF  when 
Rco  and  Rcr  are  combined  is  desired.  The  simplest  way  to  combine  them 
is  to  form  the  average  value: 
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The  PDF  for  the  sum  of  two  independent  random  variables  whose  PDFs  are 

known  is  given  by  the  convolution  integral. 7  Hence,  the  PDF  for  R^  is 

given  by  the  convolution  of  Equation  3  with  the  similar  PDF  for  Rcr 

(where  we  assume  that  y  =  y  -  y): 

co  cr 


mA>  - 


4y 


( 


(1+Y2) 


dx 


r  t~\  .  2x.  2  ,  2 

[  (1  4-  — )  +  y  (1 


2x.  2  3/2 
D  ’  1 


{[l  +|(x-ra)]2  +  y2[i  -  |(x-ra)12}3/2  j 

The  integral  in  Equation  5  does  not  have  a  closed-form  solution,  and 
thus  numerical  evaluation  is  required.  We  emphasize  that  this  form  for 
the  PDF  assumes  the  random  variables  Rcq  and  Rcr  are  independent;  if 
they  are  not,  then  a  joint  PDF  containing  the  dependence  must  be  used. 


One  way  to  arrive  at  a  PDF  for  R^  that  can  be  expressed  in  closed 
form  is  to  assume  that  the  ratio  Z  in  Equation  2  is  constant  and  that 
only  ip  is  a  random  variable.  This  simplified  case  serves  as  a  conven¬ 
ient  check  on  the  more  complicated  calculations  given  below  and  is  also 
representative  of  some  physical  situations  (see  Reference  6,  for 
example).  In  Appendix  A  we  derive  the  PDF  for  Ra  f°r  simple  averaging 
with  Z  as  a  parameter  rather  than  a  random  variable. 

An  important  point  is  the  manner  in  which  the  random  variables  in 
Equations  4  or  5  can  be  expected  to  be  related  (i.e.,  correlated)  for 
real  targets.  For  each  polarization,  the  independence  between  the 
scattered  amplitudes  and  the  phase  implicitly  assumed  for  Equation  3 
seems  to  be  firmly  established.  The  nature  of  the  correlation  between 
the  phase  of  the  co-  and  cross-polarized  echoes  for  real  targets  is 
largely  unknown;  the  same  is  true  for  the  echo  amplitude.  We  make  the 
assumption  that  the  phases  and  the  amplitudes  in  the  two  echoes  can 
reasonably  be  expected  to  have  a  correlation  ranging  from  0  to  1,  and 
that  the  correlation  between  the  phases  of  the  two  echoes  is  independent 
of  the  correlation  between  the  amplitudes  of  the  two  echoes.  Note  the 
contrast  with  the  situation  for  frequency  diversity  combining:  for 
frequency  diversity,  the  phase  (that  appears  in  Equation  2)  is  usually 
a  random  variable  from  one  frequency  to  the  next;  however,  the  amplitude 
ratio  Z  is  not,  because  the  amplitude  varies  slowly  over  the  normal 
frequency  range  used  in  frequency  diversity  combining. 


7B.  V.  Gnedenko.  The  Theory  of  Probability.  New  York,  Chelsea, 
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In  addition  to  considering  a  simple  average  of  Rco  and  Rcr,  we 
also  consider  an  average  in  which  the  values  of  Rco  and  RCf  are  weighted 
according  to  the  fractional  proportion  of  the  RCS  in  each  polarization: 


R  = 
w 


P  P 

_  C°  D  .  Cr 

P  +  P  CO  P  +  P 
co  cr  co  cr 


(6) 


where  Pco  and  Pcr  are  estimates  of  the  echo  power  in  each  polarization. 
The  rationale  for  investigating  this  combining  method  can  be  understood 
by  the  following  reasoning.  The  total  scattered  field  from  the  two 
scatterers  in,  say,  the  co-polarized  return  is  given  by 


T  2  1/2 

E  =  E„  (1  +  2Z  cosij,  +  Z  )  ' 
co  2, co  co  co  co 


(7) 


The  form  on  the  right-hand  side  of  Equation  7  also  appears  in 
Equation  2;  clearly  when  (eJo/e2,co)  is  large,  the  value  of  Rco  is 
small,  so  that  more  weight  should  be  given  to  the  Rco  or  Rcr  with  the 
higher  value  of  received  power.  Figure  3  is  a  plot  of  Equation  7  to 
the  same  scale  as  Figure  2  and  also  demonstrates  the  advantage  of  the 
weighting.  The  minimum  in  the  surface  (corresponding  to  the  lowest 
total  received  echo  power)  occurs  near  \p  =  180  degrees  and  Z  =  1.0, 
which  is  the  point  of  highest  glint  error  in  Figure  2. 


FIGURE  3.  Plot  of  Equation  7. 
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The  form  of  Equation  6  assumes  that  on  the  average  the  power  in 
the  co-  and  cross-polarized  echo  is  the  same.  The  few  experimental 
polarization  measurements8  on  aircraft  typically  indicate  that  in  fact 
the  cross-polarized  echo  power  is  about  10  to  15  dB  less  than  the  co¬ 
polarized  echo  power.  One  possible  way  to  account  for  this  difference 
in  practice  is  to  bias  the  estimate  of  Pcr  according  to 

P  AP 

Rw  =  P  +  AP  Rco  +  P  +  AP  Rcr  ’  (8 

co  cr  co  cr 


where  A  is  a  factor  to  set  the  average  value  of  Pcr  eq*  to  Pco. 
However,  we  have  not  considered  this  form  of  averaging  this  report. 


Deriving  an  equation  (even  in  integral  form)  for  F  of  Rw 

appears  to  be  a  formidable  task.  In  view  of  this,  we  1  used  Monte 
Carlo  techniques  to  obtain  plots  of  the  PDF  both  for  R  and  for  R^. 


MONTE  CARLO  SIMULATION 


We  used  10,000  trials  to  generate  the  PDFs  in  the  following  manner. 
For  each  trial,  a  pseudorandom  uniform  number  generator  on  the  interval 
0  to  1.0  generated  "draws"  for  the  angles  and  ipcr.  Although 
individual  values  for  E  could  be  generated  by  drawing  from  a  Rayleigh 
distribution  and  forming  Z,  a  faster  method  uses  the  fact  that  the  PDF 
for  the  ratio  of  two  Rayleigh-distributed  variables  is  given  by4 

f  (Z)  =  - 2y2\  ■  ■  .  (9) 

(l  +  yzT 

A  draw  from  this  distribution  can  be  accomplished  by  drawing  a  number  N 
from  a  distribution  uniform  in  the  interval  {0,1}  and  using  the  trans¬ 
formation  Z  ■  (1/y)  /(N/(N  -  1))  . 

As  discussed  in  the  Introduction,  some  degree  of  correlation  can 
be  expected  between  the  co-  and  cross-polarized  echoes.  This  corre¬ 
lation  was  accounted  for  in  the  Monte  Carlo  process  by  defining  the 
correlation  coefficient  k  as  the  fraction  of  trials  for  which  Z  is  set 

equal  to  Z  and  <p  is  set  equal  to  p  .  Thus,  for  k  =  0.5,  for 
cr  co  cr 


8I.  S.  Gradshteyn  and  I.  M.  Ryzhik.  Table  of  Integrals,  Series, 
and  Products.  New  York,  Academic  Press,  1965. 
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example,  in  5000  trials  out  of  the  10,000  total  trials,  we  set  Z  =  Z 
,  .  co 

and  u  =  ili 

co  cr 


SIMULATION  RESULTS 


We  generated  PDFs  for  Rco,  R^»  and  Rw  during  each  run  for  a  range 
of  values  of  k  and  y.  Examples  of  the  simulation  results  are  given  in 
Figures  4  through  7.  In  all  of  these  plots,  we  set  the  scatterer 
separation  D  =  1.0,  so  that  R  can  be  interpreted  as  scaled  in  units  of 
scatterer  separation. 


APPARENT  RADAR  CENTER 


FIGURE  4.  Probability  Density  Functions  as  Determined 
by  Monte  Carlo  Simulation  for  k  =  0.0.  Solid  curve,  RCS 
weighted  average  (Rw);  dot-dashed  curve,  simple  average 
(R^);  dashed  curve,  single  polarization  (Rco  or  ^cr^’ 

The  arrows  indicate  the  positions  of  the  two  scatterers. 
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APPARENT  RADAR  CENTER 

FIGURE  5.  Probability  Density  Functions  as  Determined 
by  Monte  Carlo  Simulation  for  k  =  0.5.  Solid  curve,  RCS 
weighted  average  (Rw) ;  dot-dashed  curve,  simple  average 
(Ra);  dashed  curve,  single  polarization  (Rco  or  Rcr) • 

The  arrows  indicate  the  positions  of  the  two  scatterers. 
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FIGURE  6.  Fraction  (F)  of  Area  Under  PDF  Curve  Exceeding 
Scatterer  Separation  Versus  Correlation  k. 


'.S' 


4 


FIGURE  7.  Probability  Density  Functions  as  Determined 
by  Monte  Carlo  Simulation  for  k  =  0.0  and  y  =  1.4. 
Solid  curve,  RCS  weighted  average  (Rw) ;  dot-dashed 
curve,  simple  average  (RA) ;  dashed  curve,  single 
polarization  (Rco  or  Rcr).  The  arrows  indicate  the 
position  of  the  two  scatterers. 
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Figure  4  shows  the  PDFs  for  k  =  0.0  and  y  =  1.0.  The  PDFs  for  Rw 
and  RA  both  have  lower  dispersion  than  the  PDF  for  Rco  (or  Rcr)  alone. 
The  density  is  higher  around  the  geometric  center  of  the  target  for  Rw 
and  Ra,  and  the  tails  are  lower.  For  comparing  glint  PDFs,  the  natural 
measure  of  dispersion  is  the  fraction  of  the  area  under  the  PDF  curve 
that  lies  outside  the  bounds  of  the  target,  i.e.,  that  part  of  the  PDF 
curve  satisfying  the  condition  [ R ]  >  0.5.  Denoting  this  fraction  by  F, 
the  distributions  in  Figure  4  yield  values  for  F  of  0.26,  0.18,  and  0.11 
for  Rco,  Ra,  and  Rw,  respectively.  The  combining  of  the  co-  and  cross¬ 
polarization  returns  produces  on  average  a  smaller  glint  error,  with 
the  RCS-weighted  average  producing  a  smaller  error  than  the  simple 
arithmetic  average. 
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Figure  5  shows  the  PDFs  for  k  =  0.5  and  y  =  1.0.  In  this  case  the 
decrease  in  the  spread  of  R^^  and  RA  (as  compared  with  Rco)  is  not  as 
large  as  in  Figure  4  because  of  the  correlation.  The  PDFs  for  RA  and 
Rw  in  fact  appear  quite  similar,  although  the  tails  of  the  PDF  (not 
visible  in  Figure  5)  are  substantially  smaller  and  fall  off  much  more 
rapidly  for  Rw  as  compared  with  RA.  The  values  of  F  are  0.22  and  0.19 
for  Ra  and  Rw»  respectively.  In  Figure  6  the  values  of  F  for  RA  and 
Rw  are  plotted  as  a  function  of  the  correlation  k,  displaying  a  linear 
variation  between  k  =  0  and  k  =  1.0. 
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In  Figure  7  we  plot  the  PDFs  for  k  =  0.0  and  y  =  1.4.  The  PDFs 
are  skewed  towards  the  stronger  scattering  center,  but  the  comparative 
values  of  F  (F  =  0.11,  0.20,  and  0.26  for  Rw,  R^,  and  RCo>  respectively) 
are  still  approximately  the  same  as  for  the  case  when  y  =  1.0.  The  use 
of  RCS-weighting  of  the  tracking  errors  when  combining  clearly  produces 
a  significantly  smaller  spread  in  the  glint  distribution. 


SUMMARY  AND  CONCLUSIONS 


We  have  computed  the  probability  density  function  (PDF)  for  linear 
glint  for  the  situation  in  which  the  tracking  angles  derived  from  the 
co-  and  cross-polarized  echoes  are  combined.  We  assumed  a  target  con¬ 
sisting  of  two  scattering  centers,  each  producing  an  echo  signal  with  a 
Rayleigh  distribution,  and  investigated  two  methods  of  combining  the 
tracking  angles.  Because  of  the  complexity  of  the  equations  and  the 
desire  to  account  for  correlation  between  the  co-  and  cross-polarized 
echoes,  we  generated  the  PDFs  using  a  Monte  Carlo  simulation. 

The  results  show  that  combining  the  two  tracking  angles  in  this 
manner  can  reduce  the  fraction  of  the  glint  PDF  falling  outside  the 
bounds  of  the  target  from  0.26  (the  fraction  when  either  echo  is  used 
by  itself)  to  0.11,  if  the  co-  and  cross-polarized  echoes  are  not 
correlated.  If  the  echoes  are  correlated,  the  improvement  varies 
linearly  with  the  correlation  coefficient  between  these  two  values. 
Combining  using  weighting  of  the  derived  tracking  angle  by  the  RCS  in 
each  echo  produces  a  better  improvement  than  simple  averaging  of  the 
angles . 

A  number  of  additional  features  could  conceivably  be  incorporated 
into  the  calculation.  Examples  include  the  effects  of  a  finite  response 
time  of  the  automatic  gain  control  loop  for  each  channel,  additional 
scattering  centers,  and  the  effects  of  anisotropy  of  the  scattering 
centers.  These  features  are  not  expected  to  change  the  basic  conclusions. 
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Appendix  A 

DERIVATION  OF  A  PROBABILITY  DENSITY  FUNCTION  FOR  THE  AVERAGE  GLINT 
FOR  p  AS  A  RANDOM  VARIABLE  AND  Z  AS  A  PARAMETER 


In  this  Appendix  we  derive  the  probability  density  function  (PDF) 
assuming  that  Z  is  a  parameter  and  that  only  the  relative  phase  ip  is  a 
random  variable.  Simple  averaging  of  Rco  and  Rcr  is  used;  for  simpli¬ 
city,  Zco  =  Zcr  ;  Z  is  assumed;  and  ip  is  assumed  to  be  uniformly 
distributed  on  the  interval  0  to  it. 


For  a  uniform  distribution  of  p,  the  cumulative  probability  function 
F(R)  for  either  the  co-  or  the  cross-polarized  glint  is  obtained  by 
solving  Equation  2  for  ip  and  dividing  by  n: 


F(R) 


=  —  cos 

7T 


,[i 


-  2R(1  +  Z2) 
4ZR 


(A-l) 


where  R  represents  either  Rco  or  Rcr.  The  probability  density  function 
f(R)  then  is  found  by  differentiating  Equation  A-l: 


f(R) 


_ 1 _ 

2  ttR  [  (R  -  rx)(r2  -  R)]1/2 


(A-2) 


in  which  rj  =  (0.5)  (1  -  Z) / (1  +  Z)  and  r2  *>  (0. 5)  (1  +  Z) / (1  -  Z)  are 
the  minimum  and  maximum  values,  respectively,  that  can  be  attained  by 
R  for  a  given  value  of  Z.  A  plot  of  f (R)  for  Z  =  0.5  is  shown  in 
Figure  A-l. 

The  PDF  for  R^  is  then  given  by  the  convolution  of  f(Rco)  with 
f (Rcr)  (the  factor  of  2  needed  to  form  the  average  is  taken  care  of 
after  the  calculation  of  the  convolution  by  halving  the  interval  over 
which  it  is  defined).  The  convolution  integral  is  given  by 


I  1 1  i  \f 


fTTT«T 
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FIGURE  A-l.  Method  of  Obtaining  PDF  for  Average  Glint 

by  Forming  Convolution  Integral.  The  solid  curve  is 

the  PDF  as  given  in  Equation  A-2 ,  and  the  dotted  curve 

is  this  PDF  reflected  through  the  origin  and  translated 

a  distance  R. .  The  hatched  area  is  the  convolution. 

A 


where  Rj  =  RA  -  ri  and  R2  =  RA  ~  *2-  The  limits  on  the  integral  are 
determined  by  inspection  of  Figure  A-l.  The  convolution  is  given  by 
the  common  area  shown  in  hashed  lines  as  the  mirror  image  PDF  (shown 
with  dashed  lines)  is  moved  to  the  right.  The  integration  thus  must 
be  divided  into  two  parts  with  the  limits 

R.-r,  r  9 

f(RA)  =  /  +  /  •  (A-4) 

r,  R.-r„ 


The  integral  can  be  evaluated  in  terms  of  elliptic  integrals  by 
using  the  partial  fraction  identity 
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x(£  -  x)f(x)  £  xf(x)  (£  -  x)f(x) 


(A-5) 


to  expand  Equation  A-3  into  two  integrals  and  then  applying  Equations 
3.149.4  and  3.151.4  from  Gradshteyn  and  Ryzhik.8  The  result,  after 
tedious  algebra,  is 


-  «  - ■  -  a- 


,  ••  v  .  A 
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FIGURE  A-2.  Plot  of  Equation  A-6. 

Figure  A-3  compares  a  PDF  given  by  Equation  A-6  with  the  PDF 
generated  from  the  Monte  Carlo  code  used  to  generate  the  PDFs  in  this 
report.  The  good  agreement  between  the  two  curves  serves  as  a  verifi¬ 
cation  of  the  Monte  Carlo  calculations. 


FIGURE  A-3.  Comparison  of  PDF  Obtained  from  Monte  Carlo 
Simulation  (Solid  Curve)  and  from  Equation  A-6  (Dashed 
Curve) . 
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